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Abstract 

We calculate contributions to Bg — Bg mixing through tree-level sneutrino exchange in the 
framework of the minimal supersymmetric standard model with R-parity violation, including the 
next-to-leading-order QCD corrections. We compare our results with the updated bounds on the 
Bg — Bg mass difference reported by CDF collaborations, and present new constraints on the 
relevant combinations of parameters of the minimal supersymmetric standard model with R-parity 
violation. Our results show that upper bound on the relevant combination of couplings of Bg — Bg 
mixing is of the order 10~^. We also calculate the B^ — B^ and — mass differences, and 
show that the upper bounds on the relevant combinations of couplings are two and four orders of 
magnitude stronger than ones reported in the literatures, respectively. We also discuss the case 
of complex couplings and show that how the relevant combinations of couplings are constrained 
by the updated experiment data of Bg — Bg, B^ — Bd mixing and time-dependent CP asymmetry 
Sj/ipKsi future possible observations of Sjj^^ at LHCb, respectively. 
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I. INTRODUCTION 



Very recently, the D0 collaboration and the CDF collaboration at the Fermilab Tevatron 
reported their updated measurements of the mass difference between Bg and Bs mesons 
(Ams). The new bounds on the mass difference are l|, l2|: 



D0: 17 ps"^ < Am, < 21 ps"^ 
CDF: Am, = 17.77+|]:J|] ± 0.07 ps"^ 



(1) 



It was the first time that both the lower bound and the upper bound for the Bg — Bg mixing 
are presented. Especially the CDF result has reached an accuracy of about 1%. The new 
results are important for the precision test of the standard model (SM), especially for the 
determination of the unitary triangle. Moreover, if the SM predictions are consistent with 
the above results, these data will put severe constraints on the flavor structure of the possible 
new physics models beyond the SM. 

In the literature, there have already been many discussions about the implications of the 
new measurements. In Ref. js], Q, 5, 6, 7|, the authors carried out model-independent analysis 
of the constraints on extensions of the SM. There are also model-dependent calculations in 



some new physics models beyond the SM 



11 



12 



13 



14| . In this paper, we 



further investigate the effects of the minimal supersymmetric standard model with R-parity 



violation (MSSM-RPV) 15j on the neutral meson mixing. 

R-parity is a discrete symmetry defined by Rp = (_i)3B+l+25^ where B is the baryon 
number, L is the lepton number and S is the spin of the particle. In a supersymmetric 
extension of the SM, all the particles in the SM have Rp = 1, while all the superpartners have 



Rr, 



-1. R-parity conservation is imposed in the minimal supersymmetric standard model 



(MSSM) to keep proton stable. However, this requirement is not necessary for a fundamental 
theory, and one can always introduce lepton number or baryon number violating terms in 
the Lagrangian. For the neutral meson mixing, the relevant terms in the Lagrangian are 



RPV 



(2) 



These interaction terms can induce Bg—Bg mixing through tree level sneutrino exchange, and 
thus, probably, the relevant combination of couplings ^ A^32A^23 will severely constrained 

. Similarly, the terms shown in Eq. ([2]) also can induce 
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the — and — mixing, which were discussed in Ref. 16|] at the leading-order, 



and the constraints on the relevant combinations of couplings through comparing with data 



were given by [16j 



E^^3iA>.^3.3x 10-^ 

i 

E^^2iA:>^<4.5x 10-^ (3) 

i 

where rii = (100 GeV/m^)^, m^;.^ is the mass of the sneutrino of the i-th generation. 
However, authors of Ref. 16|] did not include the contributions of the SM, since they believed 
that both theoretical and experimental results involved considerable uncertainties at that 
time. Recently, both SM theoretical predictions and experimental results has been improved 



significantly [o]- Thus, besides investigating of B^-B^ mixing, it is also worthwhile to 
reinvestigate the constraints on the combinations of R-parity violating (RPV) couplings with 
the updated data on B^-B^ and K^-K^ mixing including the SM contributions. Moreover, 
in general, the next-to-leading-order (NLO) QCD corrections are significant, so we also 
calculate the NLO QCD effects on the above neutral meson mixing in this paper. 

In addition to the above mass differences, the R-parity violating couplings can also con- 
tribute to the CP asymmetries in the meson decay processes, when the couplings are com- 
plex. So the combinations of RPV couplings discussed above can affect observables related 
to time- dependent CP violation in processes such as Bd —>■ J/ipKg and Bg —>■ J/ip<f). In 
this paper, we also consider the constraints on the relevant combinations of MSSM-RPV 
couplings from these CP violation observables. 

We organize our paper as following. Section [III is a brief summary of the formalism for 
the calculation of the neutral meson mixing , CP asymmetry in B physics and the results 
in the SM. In Section IIIII , we calculate the contributions from the MSSM-RPV at leading- 
order and next-to-leading-order in QCD. In section HVl we present our numerical results and 
discussions. 



II. BASIC FORMALISM AND ANALYTICAL RESULTS IN THE SM 

In order to make our paper self-contained, we first illustrate the basic formalism for 
the calculation of the B^i^^yB'^^^^ and K^-K^ mass differences Am, the time-dependent CP 
asymmetry in B^ decay, Sj/^Ks ^^'i the time-dependent CP asymmetry in B^ decay, Sj/^^. 
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We start from the AF = 2 effective Hamiltonian 

TYcff = ^CiQi + h.c. 

i 

The relevant operators for our concerning are 

Q2 = q'^PLb^fPRb", 



(4) 



(5) 
(6) 
(7) 



where q = d,s for Bd or Bg meson mixing, respectively. Similar expressions for — K'^ 
mixing can be obtained by replacing b with s and setting q = d. 

With this effective Hamiltonian, the Bg — Bg mass difference can be expressed as 



Amg = 2\{Bg\nos\Bg) \ =2 



Ci {Bg \ Qi \Bg) 



(8) 



aj/^Ksi<t>)y time-dependent CP asymmetry in Bd(s) — ^ J/'ipKs^cp) decays, can be expressed as 



where S 



0'j/i>K,{<t>) = 'S'j/^ii-40)sinAmrf(5)t, (9) 
sin2/?e//, and /?e// = |arg(_Brf(s)| Hcs \Bd(s))- In the framework of SM, 



SM 
J/ipKs 



SM 

J/i'<f> 



sin2/3s 



(10) 



with 



o _ {yCKM)cd{VSKM)cb -^ 

P - arg(^-— — — — j, 

\VCKM)td\ycKM)tb 
iycKM)tsiycKM)tb. 



Ps = arg( 



( VcKM ) cs ( VcKm) cb 

In the SM, only Cq is nonzero, which has been calculated to the NLO in QCD |18l . [1 
and is given by 



(11) 



-6/23 



1 + 



An 



(12) 



where rjB = 0.55 ± 0.1 20|, J5 = 1.627 in the naive dimensional regularization scheme 
(NDR), Xt = mf/m^^ and So{xt) is the Inami-Lim function 2l|. The scale /i is usually 
taken to be ~ m;, in i? physics. 





5o(Ai) 


Biifi) 




K°{H = 2GeV) 


0.69 


1.03 


0.73 


Bdifi = nib) 


0.87 


1.16 


1.91 


Bsin = nib) 


0.86 


1.17 


1.94 



TABLE I: Non-perturbative parameters used in our numerical calculations. 
The expression for the — mass difference is a httle different, which is given by 

/^rriK = 2Re {K°\n^s\K°) = 2Re ^ Q . (13) 



The nonzero NLO Wilson coefficient in the SM is 



Gj^lTlyy 

47r2 



[Xlr]iSo{xc) + \tr]2So{xt) + 2\c\tr]3So{xc, Xt)] {af\fx)) 



"2/9 



1 + ^J3 



(14) 



where Aj = V*.Vi(i, J3 = 1.895 in the NDR scheme and the parameters rji are r/i = 1.38±C 



r]2 = 0.57±0.01, r]3 = 0.47±0.04 



20 



JO, 



18| 



22l |. respectively. 5*0 functions can be found in Ref. 
The matrix elements of the operators Qi (i=0,l,2) between two hadronic states can be 
obtained by using the vacuum insertion approximation (VIA) and partial axial current 
conservation (PCAC). We refer the readers to Ref. [23| for details. The results are 

1 



{QM) = (i?.|Qo(/i) \Bs) = ^flmBBM, 



m. 



24 4(m,(^)+m,(^))2 



m 



Bs 



5i(/i), 

52(/i), 



(15) 

(16) 

(17) 



12(ms(yu) + mbdj,))"^ ^ 

where m^^ and /^^ is the mass and decay constant of the Bg mesons, respectively. Clearly, 
the expressions for Bd and mesons can be obtained by simple substitutions. -Bj(/i) is 



the non-perturbative parameters and can be calculated with lattice method 
their values to be used in our numerical calculations in table HI 



24 



25|. We list 



III. MSSM-RPV CONTRIBUTIONS AND NLO QCD CORRECTIONS 



In the MSSM-RPV, there are tree-level contributions to the -Bs(d) — -Bs(d) and K'^ — 
mixing through sneutrino exchange. The tree diagrams contribute through the operator Qi, 
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and the corresponding Wilson coefficient is 




(18) 



where jk = 32, 31, 21 for Bg, and mesons, respectively. We assume universal sneutrino 
masses for simplicity, i.e., rriy. = rriy, i = 1,2, 3. 

At NLO in QCD, both Qi and Q2 contribute due to color exchange. To calculate the 
NLO Wilson coefficients, we match the full theory onto the effective theory at the SUSY 
scale, and then run the coefficients down to the hadronic scales using the renormalization 
group equations (RGE). In our calculations, we use dimensional regularization in = 4 — 2e 
dimensions to regulate, and use MS scheme to renormalize the ultraviolet (UV) divergences. 
We keep the heaviest quark mass mq {Q = b for B mesons and Q = s for K mesons) and 
set all other quark masses to be zero in the Wilson coefficients. However, we will keep the 
lighter quark mass rriq in the intermediate stages of our calculation in order to regulate the 
infrared (IR) divergences. We also set all external momenta to zero, since the coefficients 
should not depend on them. 



The NLO diagrams in the full theory are shown in Fig. [T] In the full theory, all the UV 
divergences should be removed by the renormalization of the quark wave functions and the 
coupling constants. The renormalized amplitude in the full theory is 




g 



FIG. 1: One-loop Feynman diagrams in the full theory. 



A 



full — — 




'5 



iQi) 



tree 



(19) 
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where Cp = 4/3, = 3 is the number of colors, and (Qi)tree the tree level matrix elements 
of the operators. 




FIG. 2: The next-to-leading-order corrections in the effective theory. 



In the effective theory, there are remaining UV divergences after taking into account the 
quark field renormalization, which must be canceled by the renormalization of the effective 
operators. Calculating the diagrams in Fig. [2] with the insertion of the operators Qi and 
Q2, we get the following amplitudes after quark field renormaliztion: 

.2 

^tree, (20) 



a, / 3 rn 
l + ^C*^ T + 2-41n^ 



2 ^ a m2 ^ 
+ — In 



An 



m„ 



n a, /'3 mj^ 1 m%\ 

(Q ) = f^ - + l-31n^--ln^ (Qi) 
47r \ e /i^ N J 



tree 



+ 



4:71 N 



^lf^ + l_31n!^l+^C,|l 



2^ 



m 



In^l +pin^ 
' Air J 



a. 



{Q 2) tree ' 

(21) 

where 1/e = 1/e — 'Je + ln47r. The renormalization constant matrix for the two operators is 



, a, 1 / 8 

^ + /^ >' 
Arc e \ 3 



(22) 



from which we obtain the anomalous dimension matrix: 



7Q 



'27T 
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(23) 



Our results of anomalous dimensions agree with those in Ref. 26|. 

Matching the results in the full theory and the effective theory, we extract the Wilson 
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coefficients: 



C,(™,) = -V^fl + f^l£j^y (24) 

mi \ Atc 1 — X I 



V 



C.(™.)= ^^£^l£l!i£, (25) 

ni}, 47r A/ 1 — x 

i " 

where x = rn?Q/m1. These coefficients satisfy the renormahzation group equations 

IgMf^)) \ , (26) 



^1^^ \C^2(yU)y ^ \C2{fi) 

from which we can solve the Wilson coefficients for arbitrary scale /i. 

AniB, Sji^Ks{(t>) defined in terms of the matrix element (-Bs(d)l '^ff/"^ \Bs(d))) which 
can be written as 



{Bs{d)\ ^eff \Bs(d)) = Asm + Arpv = Asm{^ H — -. ), (27) 

•A.SM 

where Asm and Arpv denote matrix elements of the SM and the MSSM-RPV effective 
hamiltanian, respectively. With the above matrix element, Am^ and Sj/^Ks{<t>) '^^'^ 
expressed 

Arpv 



1 



A, 



■SM 



Sj/^K.i^) = Sin (^2/5(,) + arg (^1 + ^ , (28) 

where Am^^ and (3{^s) denotes the SM contributions, respectively. In Eq. fl28l) . hadronic 
uncertainties arising from hadron decay constants cancel between Asm and Arpv-, and 
hadronic uncertainties remain only in Am^^^ and As for — mixing, Auik can 
be obtained straightforwardly. 



IV. NUMERICAL RESULTS 



In this section we present our numerical results. The SUSY scale is taken to be the 
mass of the sneutrino my. The CKM matrix elements are parametrized in the Wolfenstein 
convention with four parameters A = 0.809, A = 0.2272, p = 0.197 and rj = 0.339. The other 
standard model parameters are taken to be Gp = 1.16637"^ GeV~^, as{m,z) = 0.118, mt = 
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173 GeV [27]. The mass of B meson and K meson are m^^ = 5367.5 MeV, uibq = 5279.4 MeV 
and rriK = 497.648 MeV [27|. The time-dependent CP asymmetry in Bd decay Sj/^Ks = 
0.6871q'o32 2 7|. The recent experimental values of the mass differences of the and 
mesons mixing are Am^ = (0.5292 ± 0.0009) x lO^^ pg-i, Am^, = 0.507l[]:[j[}^ ps'^ Q, 
and ArUs is shown in Eq. ([T]). 

In our numerical calculations, we first neglect the uncertainties from the hadronic param- 
eters in the SM and assume that the predictions of the SM can reproduce the central values 
of the experimental data. Thus we demand that the RPV contributions must not exceed 
the experimental upper and lower bounds of the corresponding data. 

First, we consider the new measurement of B^ — B^ mixing, which can constrain the 
combination A32 = Yl K32'^m once the sneutrino mass is given. Fig. [3] shows the allowed 

i 

region in the A32 — mc, plane, where the gray area and the dark area is the allowed region 
extracted from the leading-order amplitude. We find that with the new data, the couplings 
are constrained to the level of 10~^. After including the NLO QCD corrections, the gray 
area is excluded and only the dark area survives. One can see that since the NLO QCD 
effects increase the resulting mass difference, the bound is more stringent (about 50% lower) 
than that from LO calculations. If setting m/> = 100 GeV, the bounds on A32 from the 
calculations of LO and NLO are 

A32 < 4.6 X lO-^(LO), 2.4 X lO-^(NLO), (29) 

respectively. 

Fig. m and Fig. [5] show the constraints on the combination A31 = ^ -^^31-^^13 from Bq — Bo 

i 

mixing and A21 = X]-^i2i-^ii2 from Kq — Kq mixing. The situation here is similar to that 

_ i 

in Bs — Bs mixing: the NLO QCD corrections increase the mass differences and thus give 
more strong constraints on the combinations of couplings of MSSM-RPV. At the reference 
point mi, = 100 GeV, the bounds are 

A31 < 2.9 X 10-^°(LO), 1.5 X 10-i°(NLO), (30) 
A21 < 8.8 X lO'^^(LO), 2.8 X lO-i^(NLO). (31) 



Above bounds are two and four orders of magnitude stronger than ones given in Ref. 16| . 
as shown in Eq. ([3]), which is due to the fact that the updated data of the measurements 
and the contributions from the SM have been considered in our calculations. 
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m„-(GeV) 



FIG. 3: The bounds for ^[^^'^2^ vary with the mass of the sneutrino. The gray area and the 

i 

dark area are the allow region extracted from the LO amplitude. After NLO corrections, the gray 
area is excluded and only the dark area survives. 

We further consider the situation that the hadronic uncertainties are involved. In this 
case, the bounds on the combinations of couplings in MSSM-RPV are looser than those in 
Eqs. (12^ . (15U]1 and The new bounds from the three observables are 

A32 < 12.392 X 10"^ (LO), 

< 6.481 X 10"* (NLO); 
A31 < 7.833 X 10"^ (LO), 

< 4.096 X 10"^ (NLO); 
A21 < 3.47 X 10"^^ (LO), 

< 1.39 X 10"^^ (NLO), (32) 

respectively. 

We also discuss the situation that Aj^ is complex, and parametrize the combinations of 
relevant couplings as Ajj = |Aij|e^*'^*^ , where ij can be 32 and 31. After assuming arbitrary 
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100 150 200 250 300 

m-(GeV) 



FIG. 4: The bounds for ^ Ksi^as vary with the mass of the sneutrino. 

i 

phases aij, we can obtain constraints on the magnitude |Ajj|. Fig. Eland Fig. [7] shows the 
the allowed range of A32 and A31 for = lOOGeV, taking into account the constraints 
from Arris and Am^. For the Bg — Bg mixing, A32 is about (9(10~^), while for the B^ — B^ 
mixing, A31 is about (9(10^^^). In both figures, there are additional areas with large |Ajj| 
besides ordinary areas with small |Ajj|, which correspond to ones where contributions from 
the MSSM-RPV are larger than those from the SM, roughly twice the SM contributions but 
with different signs. We do not discuss complex couplings in K system due to the fact that 
Eq. f|T3l) holds only under the condition that imaginary part of {K^ \ Ticfr \ is far less than 
the real part. 

The combinations of complex couplings can introduce new CP-violation origins, so we 
further investigate how the CP asymmetries in B decays constrain these combinations of 



couplings. Using the experimental data of Sj/^Ks we plot the allowed areas of the 
magnitude IA31I and the phase (T31 in FigjH], which shows that the upper bound for IA31I 
is generally about (9(10^^°), except for some special values of (T31. Those special values 
correspond to arg(l + Arpv/ Asm) = 0, and we have sin(2/3efr) = sin(2/5). In fact, another 
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100 150 200 250 300 

m-(GeV) 



FIG. 5: The bounds for ^ '^i2i-^ii2 vary with the mass of the sneutrino. 

i 

area with arg(l + Arpv / Asm) = vt — 4/? also survives from constraint of experiment data of 
Sj/iiiK^y however, this possibihty has been excluded by an angular analysis of J /ipK*^ 



and a time-dependent Dalitz plot analysis of D^h^{h^ = TT^,ri^,u) 28|. For the CP 

asymmetry in Bs decays, although there is no experimental data of Sj/^^f) currently, future 
LHCb experiment will provide enough data to reach CTstatiSj/^^i,) ~ 0.03 in the first year |29| . 
which would provide a strong constraint on new physics. So we assume Sjj^^ = 0.04 ± 0.03 
and plot constraints on IA32I and (T32 in FiglHl which shows that the upper bound for IA32I 
is about O{10~^) except for some special values of (732 similar to the case of Sj/^k^- 

In conclusion, using the updated data, we have calculated contributions to Bs — Bg 
mixing, B^—Bd mixing and K^—K^ mixing in the framework of the minimal supersymmetric 
standard model with R-parity violation including NLO QCD corrections, and presented new 
constraints on the relevant couplings of MSSM-RPV. Our results show that upper bound 
on the relevant combination of couplings of B^ — Bg mixing is of the order 10~^, and upper 
bounds on the relevant combinations of couplings of Bd — Bd mixing and — mixing are 
two and four orders of magnitude stronger than ones reported in Ref. [l^, respectively. We 
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FIG. 6: The allowed range of the magnitude and phase of -^i32-^i23 constrained by the updated 

i 

experimental data of Am^. The gray area accounts for the tree level results. When the QCD 
corrections are added, the allowed area is reduced to the dark area. 

also discussed the case of complex couplings and showed that how the relevant combinations 
of couplings are constrained by the updated experiment data of Bg — Bs, B^ — Bd mixing and 
time-dependent CP asymmetry Sj/^Ks, and future possible observations of Sj/^^f, at LHCb, 
respectively. 



Note added. While preparing this manuscript the paper of 30| appeared where the 



same coupling combination from Bg — Bg mixing is also discussed. However, the authors 



of 301] mainly dealt with contributions through the box diagram. Our results induced by 



the tree-level diagram is a few orders of magnitude stronger than theirs. 
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FIG. 7: The allowed range of the magnitude and phase of SA^3^A^*3 constrained by the experi- 

i 
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